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Abstract. The role of magnetic helicity in astrophysical large-scale dynamos is 
reviewed and compared with cases where there is no energy supply and an initial 
magnetic field can only decay. In both cases magnetic energy tends to get redistributed 
to larger scales. Depending on the efficiency of magnetic helicity fluxes, the decay of 
a helical field can speed up. Likewise, the saturation of a helical dynamo can speed 
up through magnetic helicity fluxes. The astrophysical importance of these processes 
is reviewed in the context of the solar dynamo and an estimated upper limit for the 



O ' magnetic helicity flux of 10'*^ Mx^ /cycle is given. 



1. Introduction 

Self-excited dynamo action refers to the instability of a plasma in a non-magnetic 
^ I equilibrium state to amplify magnetic fields within the framework of resistive 
^ I magnetohydrodynamics (MHD). Charge separation effects are assumed absent, i.e. no 
O ' battery-type effects are explicitly involved, although they do play a role in producing a 



weak initial seed magnetic field that is needed to provide a perturbation to the otherwise 
field-free initial state. A dynamo instability may occur when the magnetic Reynolds 
number is large enough, i.e. the fluid motions and the scale of the domain are large 
' enough. This instability is normally a linear one, but some dynamos are subcritical 
and require then a finite-amplitude initial field. In the linear case one speaks about 
slow and fast dynamos depending on whether or not the growth rate of the dynamo 
scales with resistivity. For fast dynamos the growth rate scales with the rms velocity 
of the flow, which is turbulent in most cases. Dynamos saturate when the magnetic 
energy becomes comparable with the kinetic energy. Some of the kinetic energy is then 
channelled through the magnetic energy reservoir and is eventually dissipated via Joule 
heating. 

There are several applications where one considers non self-excited dynamo action. 
Examples can be found in magnetospheric physics and in plasma physics where one is 
interested in the electromotive force induced by a flow passing through a given magnetic 
field. Later in this paper we will discuss the reversed field pinch (RFP) because of its 
connection with the a effect that plays an important role in astrophysical dynamos. The 



Magnetic helicity in astrophysical dynamos 



2 



physics of the RFP has been reviewed by Ortolani & Schnack (1993) and, in a broader 
context with reference to astrophysical large-scale dynamos, by Ji & Prager (2002) and 
Blackman & Ji (2006). 

The a effect is one of a few known mechanisms able to generate large-scale magnetic 
fields, i.e. fields whose typical length scale is larger than the scale of the energy carrying 
motions. It was also one of the first discussed mechanisms able to produce self-excited 
dynamo action at all. Indeed, Parker (1955) showed that the swirl of a convecting flow 
under the influence of the Coriolis force can be responsible for producing a systematically 
oriented poloidal magnetic field from a toroidal field. The toroidal field in turn is 
produced by the shear from the differential rotation acting on the poloidal field. Parker's 
paper was before Herzenberg (1958) produced the first existence proof of dynamos. Until 
that time there was a serious worry that Cowling's (1933) anti-dynamo theorem might 
carry over from two-dimensional fields to three-dimensional fields, as is evident from 
sentimental remarks made by Larmor (1934). 

Larmor (1919) proposed the idea of dynamo action in the astrophysical context 
nearly 100 years ago. Nowadays, with the help of computers, it is quite easy 
to solve the induction equation in three dimensions in simple geometries and 
obtain self-excited dynamo solutions with as little as 16^ mesh points using, 
for example, the sample "kin-dynamo" that comes with the PENCIL CODE 
(http : //pencil-code . googlecode . com). 

In addition to dynamos in helical flows, which can generate large-scale fields, 
there are also dynamos in non-helical flows that produce only small-scale fields. This 
possibility was first addressed by Batchelor (1950) based on the analogy between the 
induction equation and the vorticity equation. Again, this was not yet very convincing 
at the time. The now accepted theory for small-scale dynamos was first proposed by 
Kazantsev (1968) and the first simulations were produced by Meneguzzi et al. (1981). 
Such simulations are computationally somewhat more demanding and require at least 
64^ mesh points (or collocation points in spectral schemes). In the past few years this 
work has intensified (Cho et al. 2002, Schekochihin et al. 2002, Haugen 2003). We 
will not discuss these dynamos in the rest of this paper. Instead, we will focus on 
large-scale dynamos. More specifically, we focus here on a special class of large-scale 
dynamos, namely those where kinetic helicity plays a decisive role (the so-called a 
effect dynamos). Nevertheless, we mention at this point two other mechanisms that 
could produce large-scale magnetic field without net helicity. One is the incoherent a- 
shear effect that was originally proposed by Vishniac & Brandenburg (1997) to explain 
the occurrence of large-scale magnetic fields in accretion discs, and later also for other 
astrophysical applications (e.g.. Proctor 2007). It requires the presence of shear, because 
otherwise only small-scale magnetic fields would be generated (Kraichnan 1976, Moffatt 
1978). The other mechanism is the shear-current effect of Rogachevskii & Kleeorin 
(2003), which can operate if the turbulent magnetic diffusion tensor is anisotropic, so 
the mean electromotive force from the turbulence is given by —TjijJj such that the sign 
of rjijUij is positive, and that this quantity is big enough to overcome resistive effects. 
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Here, a comma denotes partial differentiation, U is the mean flow, J is the mean current 
density, and summation over repeated indices is assumed. This effect too requires shear, 
because otherwise rjijUij would be zero. Simulations show large-scale dynamo action 
in the presence of just turbulence and shear, and without net helicity, but there are 
indications that this process may also be just the result of incoherent a-shear dynamo 
action (Brandenburg et al. 2008a). 

Throughout this paper, overbars denote suitable spatial averages over one or two 
coordinate directions. Furthermore, we always assume that the scale of the energy- 
carrying eddies is at least three times smaller than the scale of the domain. We refer 
to this property as "scale separation". In this sense, scale separation is a natural 
requirement, because we want to explain the occurrence of fields on scales large compared 
with the scale of the energy- carrying motions. Scale separation has therefore nothing 
to do with a gap in the kinetic energy spectrum, as is sometimes suggested. 

Much of the work on a effect dynamos has been done in the framework of analytic 
approximations. However, this is only a technical aspect that is unimportant for 
the actual occurrence of large-scale fields under suitable conditions. This has been 
demonstrated by numerical simulations, as will be discussed below. 

2. Helical leirge-scale dynamos 

A possible way of motivating the physics behind helical dynamo action is the relation 
to the concept of an inverse turbulent cascade. This particular idea was first proposed 
by Prisch et al. (1975) and is based on the conservation of magnetic helicity, 

Hm^ f A-BdV, (1) 

J V 

where A is the magnetic vector potential and B — 'V x A is the magnetic field in a 
volume V. 

It is convenient to define spectra of magnetic energy and magnetic helicity, Euik) 
and Huik), respectively. As usual, these spectra are obtained by calculating the three- 
dimensional Fourier transforms of magnetic vector potential and magnetic field, A^ 
and Bk, respectively, and integrating |.Bfcp and the real part of A^. ■ B^. over shells of 
constant A; = |fc| to obtained Euik) and Huik), respectively. (Here, an asterisk denotes 
complex conjugation.) These spectra are normalized such that / EM{k)dk — (B^)/2/io 
and / Hyiik) dk — {A • B) for k from to oo, where angular brackets denote volume 
averages over a periodic domain and /xq is the vacuum permeability. Using the Schwartz 
inequality one can then derive the so-called realizability condition, 

k\HM{k)\/2n^<EM{k). (2) 

For fully helical magnetic fields with (say) positive helicity, i.e. Hyi — 2^QEy[{k) /k, one 
can show that energy and magnetic helicity cannot cascade directly, i.e. the interaction 
of modes with wavenumbers p and q can only produce fields whose wavevector k — p+q 
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Wave number 

Figure 1. Magnetic energy spectra at different times (increasing roughly by a factor 
of 2). The curve with the right-most location of the peak corresponds to the initial 
time, while the other lines refer to later times (increasing from right to left). Note 
the propagation of spectral energy to successively smaller wavenumbers fc, i.e. to 
successively larger scales. Adapted from Christensson et al. (2001). 



has a length that is equal or smaller than the maximum of either \p\ or \q\ (Frisch et 
al. 1975), i.e. 

\k\ < max(|p|, |q|). (3) 

This means that magnetic helicity and magnetic energy are transformed to progressively 
larger length scales. A clear illustration of this can be seen in decaying helical turbulence. 
Figure [1] shows magnetic energy spectra from a simulation of Christensson et al. (2001) 
at different times for a case where the initial magnetic field was fully helical and 
had a spectrum proportional to with a resolution cutoff near the largest possible 
wavenumber. Note that the entire spectrum appears to shift to the left, i.e. toward larger 
length scales, in an approximately self-similar fashion. The details of the argument that 
led to Equation ([3]) are due to Frisch et al. (1975), and can also be found in the reviews 
by Brandenburg et al. (2002) and Brandenburg & Subramanian (2005a). 

In the non-decaying case, when the flow is driven by energy input at some forcing 
wavenumber fcf, the inverse cascade is clearly seen if there is sufficient scale separation, 
i.e. if kf is large compared with the smallest wavenumber ki that fits into a domain of 
size L = 2TT/ki. An example is shown in Figure [2], where kinetic energy is injected at 
the wavenumber kf = 30ki. It is evident that there are two local maxima of spectral 
magnetic energy, one at the forcing wavenumber fcf, and another one at a smaller 
wavenumber that we call /Cm, which is near 7ki in Figure [2j During the kinematic stage 
the entire spectrum moves upward, with the spectral energy increasing at the same rate 
at all wavenumbers. Eventually, when the field has reached a certain level, the spectrum 
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Figure 2. Magnetic energy spectra for a run with forcing at = 30. The times, 
in units of {cski)~^, range from (dotted Une) to 10, 30, 290 (sohd hues). The 
thick sohd hnc gives the final state at Cgfcii = 1000, corresponding to Uimsfcft « 2000 
turnover times. Here, Cg is the sound speed and Urms is the turbulent rms velocity. 
Note that at early times the spectra peak at fcmax ~ 7fci. The and slopes 
are given for orientation as dash-dotted hues. Adapted from Brandenburg (2001). 



begins to change its shape and the second local maximum at moves toward smaller 
values, suggestive of an inverse cascade. However, a more detailed analysis (Brandenburg 
2001) shows that the energy transfer is nonlocal, i.e. most of the energy is transferred 
directly from the forcing wavenumber to the wavenumber where most of the mean field 
resides. This suggests that we have merely a nonlocal inverse transfer rather than a 
proper inverse cascade, where the energy transfer would be local in spectral space. 

The position of the local maximum can readily be explained by mean-field dynamo 
theory with an a effect. The evolution equation of such a dynamo is 

OB 

— = VxaB + r]TV^B, (4) 

where a is a pseudo-scalar, rjT = rj + rjt is the sum of microscopic Spitzer resistivity and 
turbulent resistivitjjjl, and an overbar denotes a suitably defined spatial average (e.g. 
planar average). Assuming B = B^exp^Xt + ik ■ x) with eigenfunction B^, one finds 
the dispersion relation to be (Moffatt 1978) 

X{k) = \a\k — r]Tk'^ , (5) 



I Note that resistivity and magnetic diffusivity differ by a /io factor. Here, we always mean the magnetic 
diffusivity, although we use the two names sometimes interchangeably. 
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where k = \k\. The maximum growth rate is attained for a value of k where dX/dk = 0, 
i.e. for k = k^ = a/27]^. The migration of the spectral maximum to smaller k can 
then be explained as the result of a suppression of a. In other words, as the dynamo 
saturates, a decreases, and so does k^, which corresponds to the spectral maximum 
moving to the left. 

It should be noted that there is also the possibility of a suppression of rit, which 
would work the other way, so that this interpretation might not work. Indeed, there 
are arguments for a suppression of r]t that would be as strong as that of a, but this 
applies only to the two-dimensional case and has to do with the conservation of the 
mean-squared vector potential in that case (Gruzinov & Diamond 1994). Simulations 
also find evidence that in three dimensions the suppression of a is stronger than that of 
T]t (Brandenburg et al. 2008b). 

In the following we discuss in detail the role played by magnetic helicity. This 
has been reviewed extensively in the last few years (Ji & Prager 2002, Brandenburg & 
Subramanian 2005, Blackman & Ji 2006). 

3. Slow saturation 

In a closed or periodic domain, the saturation of a helical large-scale dynamo is found 
to be resistively slow and the final field strength is reached with a time behavior of the 
form (Brandenburg 2001) 



where B^q is the equipartition field strength and ts is the time when the slow saturation 
phase begins. We emphasize that it is the microscopic r] that enters Equation (j6]), and 
that the relevant length scale, 2n/k^, is that of the large-scale field, so the saturation 
behavior is truly very slow. 

The reason for this slow saturation behavior is related to the conservation of 
magnetic hehcity, which obeys the evolution equation (e.g., Ji et al. 1995, Ji 1999) 



where F-^ is the magnetic helicity flux, but for the periodic domain under consideration 
we have V ■ F-^ = 0. Clearly, in the final state we have then {J ■ B) = 0. This can only 
be satisfied for nontrivial helical fields if small-scale and large-scale fields have values 
of opposite sign, but equal magnitude, i.e. (j ■ b) = —{J-B), where B = B -\- b and 
J = J + j are the decompositions of magnetic field and current density into mean 
and fluctuating parts. Here we choose to define mean fields as one- or two-dimensional 
coordinate averages. Examples include planar averages such as xy, yz, or xz averages 
in a periodic Cartesian domain, as well as one-dimensional averages such as ?/ or <^ 
averages in Cartesian or spherical domains, (r, 6, </>), where the other two directions are 
non-periodic. 



{W) (t) = Bl ^ fl - e-^"^- for t > t 



(6) 



^(A-B) = -2r^fio{J-B)-{V-Fn). 



(7) 
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Equation ([2D can be derived under the assumption that large-scale and small-scale 
fields are fully helical with B = k^A ■ B = TkmB"^ and {j ■ b) = ±kf{b'^) ^ ±k{B^^, 
where upper and lower signs refer to positive and negative helicity of the small- 
scale turbulence and we have assumed that the small-scale field has already reached 
saturation, i.e. (6^) ~ B"^^ = /io(pw^), where p is the density and u = U — U is the 
fluctuating velocity. 

A mean-field theory that obeys magnetic helicity conservation was originally 
developed by Kleeorin & Ruzmaikin (1982) and has recently been applied to explaining 
slow saturation (Field & Blackman 2002, Blackman & Brandenburg 2002, Subramanian 
2002). The main idea is that the a effect has two contributions (Pouquet et al. 1976), 

a = ax + dMi (8) 

where = — |ra; • u is the usual kinetic a effect related to the kinetic helicity, with 
ci; = V X n being the vorticity, and au = ^tJ ■ b/ po is a magnetic a effect that can, for 
example, be produced by the growing magnetic field in an attempt to conserve magnetic 
helicity. (Here, po is an average density, but we note that there is at present no adequate 
theory for compressible systems with nonuniform density.) 

Note that om is related to the small-scale current helicity and hence to the 
small-scale magnetic helicity which, in turn, obeys an evolution equation similar to 
Equation ([7]), but with an additional production term, 2£ ■ B, that arises from mean- 
field theory via £ = aB — ritpoJ, and thus from the a effect itself. The equation for 
ttM has then the form 

where i?m = Vt/v is a measure of the ratio of turbulent to microscopic magnetic 
diffusivity. With rjt = Urms/3A;f (Sur et al. 2008) we can relate this to the more 
usual definition for the magnetic Reynolds number, = Mrms/''7^f; 

via — 3i?]-Q- 

Furthermore, we have allowed for the possibility of fiuxes of magnetic and current 
helicities that also lead to a fiux of om- Such fiuxes are primarily important in 
inhomogeneous domains and especially in open domains where one can have an outward 
hehcity fiux (Ji 1999). 

The use of Equation ([8]) is sometimes criticized because it is based on a closure 
assumption. Indeed, there are questions regarding the meaning of the term j ■ b and 
whether it really applies to the actual field, or the field in the unquenched case. This has 
been discussed in detail in a critical paper by Radler & Rheinhardt (2007). Part of this 
ambiguity can already be clarified in the low conductivity limit. Sur et al. (2007) have 
shown that one can express a either completely in terms of the helical properties of the 
velocity field or, alternatively, as the sum of two terms, a so-called kinetic a effect and 
an oppositely signed term proportional to the helical part of the small scale magnetic 
field. However, it is fair to say that the problem is not yet completely understood. The 
strongest argument in favor of Equations ([H]) and ([9]) is that they reproduce catastrophic 
(i.e. i?m-independent) quenching of a and that this approach has led to the prediction 
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that such quenching can be alleviated by magnetic helicity fluxes. This prediction has 
subsequently been tested successfully on various occasions (Brandenburg 2005, Kapyla 
et al. 2008). Finally, it should be noted that Equation ([8]) has also been confirmed 
directly using turbulence simulations (Brandenburg et al. 2005c, 2007). 

The idea to model dynamo saturation and suppression of a by solving a dynamical 
equation for is called dynamical quenching. In addition to the resistively slow 
saturation behavior described by Equation IQ, this approach has also been applied to 
decaying turbulence with helicity (Yousef et al. 2003; Blackman & Field 2004), where 
the conservation of magnetic helicity results in a slow-down of the decay. This can be 
modeled by an a effect that offsets the turbulent decay proportional to rjtk"^ such that 
the decay rate becomes nearly equal to the resistive value, r]k^. This is explained in 
detail in the following section. 



4. Decay in a Cartesian domain 

In the context of driven turbulence, the properties of solutions of a decaying helical 
magnetic field were studied earlier by Yousef et al. (2003), who found that for fields 
with B^/B^^ ^ the decay of B is slowed down and can quantitatively be described 
by the dynamical quenching model. This model applies even to the case where the 
turbulence is nonhelical and where there is initially no a effect in the usual sense. 
However, the magnetic contribution to a is still non- vanishing, because the om term is 
driven by the helicity of the large-scale field. 

To demonstrate this quantitatively, Yousef et al. (2003) have adopted a one-mode 
approximation with B = B{t) exp{ikiz), and used the mean-field induction equation 
together with the dynamical a-quenching formula (JH]), 

^ = ik,x£~7^klB, (10) 
at 

where the fiux term is neglected, £ = aB — rftiki x B is the electromotive force, and 
fci = (0,0,fci). _ 

Figure [3] compares the evolution of B/B^q for helical and nonhelical initial 
conditions, B oc (l,i, 0) and B oc (1,0,0), respectively. In the case of a nonhehcal 
field, the decay rate is not quenched at all, but in the helical case quenching sets in 
for B'^/B^^ ^ R^. The onset of quenching at B"^ / Bl^ ^ R^^^ is well reproduced by 
the simulation. In the nonhelical case, however, some weaker form of quenching sets in 
when B"^ / Bl^ ^ 1. We refer to this as standard quenching (e.g. Kitchatinov et al. 1994) 
which is known to be always present; see Equation ([6]). Blackman & Brandenburg (2002) 
found that, for a range of different values of i?m, g = ?> results in a good description of 
the simulations of cyclic afi-type dynamos that were reported by Brandenburg et al. 
(2002.). 
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Figure 3. Dynamical quenching model with helical and nonhelical initial fields, 
and an additional r]t quenching fimction, rit — fyto/(l + .9l-B|/i3oq)- The quenching 
parameters are 5 = (solid line) and 3 (dotted line). The graph for the nonhelical 
cases has been shifted in t so that one sees that the decay rates are asymptotically 
equal at late times. The value of rj^ used to normalize the abscissa is based on the 
unquenched value. Adapted from Yousef et al. (2003). 



5. Relevance to the reversed field pinch 

The dynamical quenching approach has been applied to modeling the dynamics of the 
reversed field pinch, where one has an initially helical magnetic field of the form 



where we have adopted cylindrical coordinates {r,9,z). In a cylinder of radius R such 
an initial field becomes kink unstable when kR ^ vr. Both laboratory measurements 
(e.g. Caramana & Baker 1984) and numerical simulations (Ho et al. 1989) confirm the 
idea that the field-aligned current leads to kink instability, and hence to small-scale 
turbulence and thereby to the emergence of % and a. 

Just as in the case discussed in Section HI the emergence of a slows down the decay 
in such a way that the toroidal field remains nearly constant and is maintained against 
resistive decay; see Figure HI The details of this mechanism have been discussed by 
Ji & Prager (2002). In particular, they show that the parallel electric field cannot be 
balanced by the resistive term alone, and that there must be an additional component 
resulting from small-scale correlations of velocity and magnetic field, u x b, that explain 
the observed profiles of mean electric field and mean current density. Furthermore, the 
parallel electric field reverses sign near the edge of the device, while the parallel current 
density does not. Again, this can only be explained by additional contributions from 
small-scale correlations of velocity and magnetic field. The RFP experiment also shows 
that magnetic helicity evolves on time scales faster than the resistive scale, which is only 
compatible with the presence of a finite magnetic helicity fiux divergence (Ji et al. 1995, 




(12) 



Ji 1999). 
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Figure 4. Time evolution of toroidal flux in the RFP experiment of Caramana & 
Baker (1984) compared with a calculation with no dynamo effect. The decay around 
t = 18ms is due to termination of the applied electric field. Courtesy of E. J. Caramana 
and D. A. Baker. 

6. Magnetic helicity in realistic dynamos 

Astrophysical dynamos saturate and evolve on dynamical time scales and are thus not 
resistively slow. Current research shows that this can be achieved by expelling magnetic 
helicity from the domain through helicity fluxes. This is why we have allowed for the 
V ■ Fa term in Equation ([9]). Since the magnetic a effect is proportional to the current 
helicity of the fluctuating fleld, the flux should be proportional to the current helicity 
flux of the fluctuating fleld. 

The presence of the flux term generally lowers the value of | | , and since the a m 
term quenches the total value of a (= + om), the effect of this helicity flux is to 
alleviate an otherwise catastrophic quenching. Indeed, in an open domain and without 
a flux divergence, the au/ Rm term in Equation ([9]) can result in "catastrophically low" 
saturation fleld strengths that are by a factor RH"^ smaller than the equipartition fleld 
strength (Gruzinov & Diamond 1994; Brandenburg & Subramanian 2005b). 

Magnetic helicity obeys a conservation law and is therefore conceptually easier to 
tackle than current helicity. However, there is the difficulty of gauge dependence of 
magnetic helicity density and its flux. It is therefore safer to work with the current 
helicity, which is also the quantity that enters in Equations (jHl) and ([9]). However, more 
work needs to be done to establish the connection between the two approaches. 

Over the past 10 years there has been mounting evidence that the Sun sheds 
magnetic helicity (and hence current helicity) through coronal mass ejections and other 
events. Understanding the functional form of such fluxes is very much a matter of 
ongoing research (Subramanian & Brandenburg 2004, 2006, Brandenburg et al. 2009). 
In the following section we present a simple calculation that allows us to estimate the 
amount of magnetic helicity losses required for the solar dynamo to work. 
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7. Estimating the required magnetic helicity losses 

In order to estimate the magnetic helicity losses required to alleviate catastrophic 
quenching we make use of the relation between the a effect and the divergence of the 
current helicity flux (Brandenburg & Subramanian 2005a), 

" = 1 + fl„SVB?, • 

where a = da / dt and Fq is the mean flux of current helicity from the small-scale field, 
(V X e) X (V X 6), where e is the fiuctuating component of the electric field; see also 
Subramanian & Brandenburg (2004). In the steady-state limit and at large Rm we have 

"^''^ W 2kjW ^^^^ 

We neglect the J ■ B term, because the catastrophic quenching in dynamos with 
boundaries has never been seen to be alleviated by this term (Brandenburg & 
Subramanian 2005b). Thus, we use Equation f|T^ to estimate V ■ Fq as 

V ■ Fc = 2akJW. (15) 

Next, we take the volume integral over one hemisphere, i.e. 

Cc=i Fc-dS= fv-FcdV= '^R'{2aklW), (16) 

where Cq is the "luminosity" or "power" of current helicity. We estimate (a) using afl 
dynamo theory which predicts that (e.g., Robinson & Durney 1982, see also Brandenburg 
& Subramanian 2005a) 

ak.An ^ uj%^, (17) 

where uj^jc = 27r/Tcyc is the cycle frequency of the dynamo, Tcyc is the 22 year cycle 
period of the Sun, a is assumed constant over each hemisphere, and AO is the total 
latitudinal shear, i.e. about 0.3f2 for the Sun. There is obviously an uncertainty in 
relating local values of a to volume averages. However, if we do set the two equal, we 
obtain at least an upper limit for Cq. We may then relate this to the luminosity of 
magnetic helicity, £h, that we assume to be proportional to £c via 

Cc = kjC^. (18) 

With this we find for the total magnetic helicity loss over half a cycle (one 11 year cycle) 

where we have estimated ki = 27r/L for the relevant wavenumber of the dynamo in 
terms of the thickness of the convection zone L. Inserting now values relevant for the 
Sun, L = 200 Mm, R = 700 Mm, uJcyc/AVl = lO'^, and B ~ 300 G, we obtain 

i^nTcyc < 10^^ MxVcycle. (20) 
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This is comparable to earlier estimates based partly on observations (Berger & 
Ruzmaikin 2000; DeVore 2000) and partly on turbulence simulations (Brandenburg 
& Sandin 2004), but we recall that Equation ( l20i) is only an upper limit. Furthermore, 
the connection between current helicity and magnetic helicity assumed in relation (ITSll 
is quite rough and has only been seriously confirmed under isotropic conditions. In 
addition, it is not clear that the gauge-invariant magnetic helicity flux defined by 
Berger & Field (1984) and Finn & Antonsen (1985) is actually the quantity of interest. 
Following earlier work of Subramanian & Brandenburg (2006), the gauge- invariant 
magnetic helicity is not in any obvious way related to the current helicity which, in 
turn, is related to the density of the flux linking number and is at least approximately 
equal to the magnetic helicity in the Coulomb gauge. 

In summary, although magnetic helicity is conceptually advantageous in that it 
obeys a conservation equation, the difficulty in dealing with a gauge- dependent quantity 
can be quite serious. Moreover, as emphasized before, it is really the current helicity that 
is primarily of interest, so it would be useful to shift attention from magnetic helicity 
fluxes to current helicity fluxes. 

8. Conclusions 

In this review we have attempted to highlight the importance of magnetic helicity in 
modern nonlinear dynamo theory. Much of the early work on the reversed field pinch 
since the mid 1980s now proves to be extremely relevant in view of the possibility of 
resistively slow saturation by a self-inflicted build-up of small-scale current helicity, j ■ b, 
as the dynamo produces large-scale current helicity, J ■ B. Since this concept is not yet 
universally accepted, the additional evidence from the reversed field pinch experiment 
can be quite useful. 

The interpretation of resistively slow saturation has led to the proposed solution 
that magnetic helicity fluxes (or current helicity fluxes) are responsible for removing 
excess small-scale magnetic helicity from the system. This allows the dynamo to reach 
saturation levels that can otherwise be i?^^ times smaller than the equipartition value 
given by the kinetic energy density of the turbulent motions. The consequences of this 
prediction have been tested in direct simulations (Fig. 3 of Brandenburg 2005 and Fig. 17 
of Kapyla et al. 2008), confirming thereby ultimately basic aspects of incorporating the 
magnetic helicity equation into mean-field models. However, more progress is needed in 
addressing questions regarding the relative importance of magnetic and current helicity 
fluxes through the surface compared to diffusive fluxes across the equator (Brandenburg 
et al. 2009). 

We reiterate that the reversed field pinch experiment is not directly relevant to the 
self-excited dynamo, but rather its nonlinear saturation mechanism. So far, successful 
self-excited dynamo experiments have only been performed with liquid sodium (Gailitis 
et al. 2000; Stieglitz & Miiller 2001; Monchaux et al. 2007). This may change in future 
given that one can usually achieve much higher magnetic Reynolds numbers in plasmas 
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than in liquid metals (Spence et al. 2009). It might then, for the first time, be possible 
to address experimentally questions regarding the relative importance of magnetic and 
current helicity fluxes for dynamos. 
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